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Abstract 

A unified approach has been developed to study nonlinear dynam- 
ics of a ID lattice of particles with long-range power-law interaction. 
A classical case is treated in the framework of the generalization of 
the well-known Frenkel-Kontorova chain model for the non-nearest 
interactions. Quantum dynamics is considered following Davydov's 
approach for molecular excitons. 

In the continuum limit the problem is reduced to dynamical equa- 
tions with fractional derivatives resulting from the fractional power 
of the long-range interaction. Fractional generalizations of the sine- 
Gordon, nonlinear Schrodinger, and Hilbert-Schrodinger equations have 
been found. 

There exists a critical value of the power s of the long-range po- 
tential. Below the critical value (s < 3, s ^ 1, 2) we obtain equations 
with fractional derivatives while for s > 3 we have the well-known 
nonlinear dynamical equations with space derivatives of integer order. 

Long-range interaction impact on the quantum lattice propagator 
has been studied. We have shown that the quantum exciton propaga- 
tor exhibits transition from the well-known Gaussian-like behavior to 
a power-law decay due to the long-range interaction. A link between 
ID quantum lattice dynamics in the imaginary time domain and a 
random walk model has been discussed. 

PACS number(s): 03.65.-w, 03.65. Db, 05.30.-d, 05.40. Fb 
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1 Introduction 

Discrete nonlinear lattice models have become one of the most popular mod- 
els of nonlinear physics, often used to investigate a rather broad set of physical 
phenomena and systems. The first simplest ID lattice model was introduced 
a long time ago by Frenkel and Kontorova |T to study the structure and dy- 
namics of dislocations in metals (for more details see review j2] and references 
therein). The Frenkel and Kontorova discrete classical dynamics is described 
as a ID chain of atoms with the nearest-neighbor intersite/interatomic in- 
teraction and a periodic on-site potential. Despite its simplicity the model 
has many applications in solid-state and nonlinear physics, including propa- 
gation of charge-density waves, the dynamics of absorbed layers of atoms on 
crystal surfaces, commensurable-incommensurable phase transitions, and do- 
main walls in magnetically ordered structures, see for details Moreover 
the model is surprisingly attractive because it provides exactly integrable 
cases in discrete and continuous approximations including the well-known 
sine-Gordon nonlinear equation (see, for example (3]). 

A specific, but not less important application, is the case of long-range 
interatomic interaction, which is the main subject of this paper and which has 
been studied for different physical systems. Examples include ferromagnetic 
chains |3], exciton transfer in molecular crystals (§], commensurable- 
incommensurable phase transitions [7j, theory of Josephson junctions jH], [§]■ 
To model the long-range interaction potential in ID lattice, few potentials 
have been applied: power-law interaction [TOj, Lennard- Jones long-range 
coupling ^T] and exponential interaction of Kac-Beker form ^2] ■ 

Since the Fermi, Pasta and Ulam work analytical developments and 
numerical simulation on discrete nonlinear dynamical equations have pro- 
vided numerous types of nonlinear excitations with various properties: soli- 
tons, continuous and discrete breathers [2], [Ej, [IS], self-trapping states [Trj] 
and others. 

Another broad area of applications is related to the coupled system of os- 
cillators with long-range interactions and the phenomenon of synchronization 
in the system [T7]-[T§]. In many situations, impact of long-range coupling 
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can be compared to a phase transition [2Dl-j22]- Long-range coupling can be 
presented as nonlocality with finite space scale a, when the coupling energy 
is proportional to exp(— \x n — x m \/a) for atoms located at the points x n and 
x m , or with scale free interaction when the coupling energy is proportional to 
\x n — Xm\~ s (s > 0). In addition to the well known interaction with integer 
values of s , some complex media can be described by fractional values of s 
(see, for example, references in 

The goal of this paper is to focus on analytical developments on classical 
and quantum ID lattice dynamics with fractional power-law long-range in- 
teraction. The long-range interaction leads in general to a nonlocal integral 
term in the equation of motion if we go from discrete to continues space. We 
show that in some particular cases for power-law interaction with non-integer 
power s the integral term can be expressed through a fractional order deriva- 
tive. In other words, nonlocality originating from the long-range interaction 
is revealed in the dynamics in the form of space derivatives of fractional or- 
der. The appearance of fractional differential equations in the continuum 
space limit of lattice dynamics leads to a possibility to apply powerful tools 
of random processes theory [23], |2S] and fractional calculus [2E]-[2S| to the 
lattice dynamics and kinetics. Particularly, it helps to obtain some quali- 
tative results based on the known features of fractional derivatives. As an 
example, let us mention that the fractional oscillator appears to have power 
law dissipation at t — > oo and, being perturbed, gives rise to a new type of 
stochastic attractor (231- I n the case of a fractional space derivative, one can 
expect new types of synchronizations and coherent structures j3U|; where our 
approach to obtain equations with fractional derivative is applied to discrete 
model of oscillating media with long-range power-type coupling between os- 
cillators. It also worthwhile to mention that fractional dynamical equations 
prove to be an adequate approach to study chaotic dynamical systems and 
their anomalous transport properties |3*Tj . 

We present new developments on a ID lattice model with fractional 
power-law interatomic interaction defined by fractional values of the pa- 
rameter s and non-linear self-interaction potential for classical and quan- 
tum mechanical consideration. We show that the dynamics of the ID lat- 
tice can be equivalently presented by the corresponding fractional nonlinear 
equation in the long-wave limit. We concentrate on the conditions of such 
equivalence, type of the equations with fractional derivatives and some re- 
lated properties. As example of our developments fractional sine-Gordon and 
wave-Hilbert nonlinear equations have been found for classical lattice dynam- 
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ics, and fractional nonlinear Schrodinger and nonlinear Hilbert-Schrodinger 
equations have been obtained for quantum lattice dynamics in the long wave 
limit. 

The paper is organized as follows. In Sec. 2 we study classical ID lattice 
non-linear dynamics. It has been found that depending on the long-range 
interaction parameter s, the interaction term can be presented in long-wave 
limit in differential or fractional differential form which results either as a 
nonlinear wave equation or as a nonlinear fractional wave equation. Examples 
are fractional sine-Gordon and wave-Hilbert nonlinear equations. 

Sec. 3 focuses on quantum ID lattice dynamics. Depending on the long- 
range interaction parameter s, the interaction term can be presented in the 
long-wave limit in differential or fractional differential form which results 
in the nonlinear Schrodinger equation or in a nonlinear Hilbert-Schrodinger 
equation. 

Quantum lattice propagator is introduced in Sec. 4. It has been observed 
that depending on the parameter s the quantum lattice propagator exhibits 
a transition from the well-known Gaussian-like behavior to power-law decay 
due to long-range interaction. We show a link between the quantum lattice 
dynamics and random walk in the imaginary time domain. 

In conclusion we outline the new developments and discuss their applica- 
bility to other fields. 

In Appendices we briefly review the definition of the polylogarithm and 
its integral and power series representations. 



2 Classical nonlinear lattice dynamics 

Nonlinear classical discrete dynamics of a ID infinite lattice can be described 
by Hamiltonian function H(u), which depends on particle ("atom") displace- 
ment u n {t) at the site n 



H{u) = -M £ u 2 n + - £ £ V„{u n - u m f + £ U(u n ), (1) 

n=~ oo n=-oom=-oo n=— oo 

n^ra 

where the first term is the kinetic energy of the chain, M is particle mass, and 
V n - m is interaction matrix which describes long-range interaction between 
sites n and m 
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Vn—m — i r ■ (2) 

\n — m\ s 

Here Vq is interaction constant, parameter s covers different physical mod- 
els, and U(u) characterizes nonlinear interaction of the lattice atoms with 
external on-site potential. Integer values of s can be used to describe the 
well-known physical models: the nearest-neighbor approximation (s = oo), 
the dipole-dipole interaction (s = 3), the Coulomb potential (s = 1). 

Our main interest will be in fractional values of s that can appear for 
more sophisticated interaction potentials attributed to complex media. The 
equation of motion follows from Eq.([T|) 

00 dU(u) 

Mu n + Yl Vn-miUn - U m ) + = 0, (3) 

m=-oo 

To go from the discrete version Eq.(JHJ) to the continuum one, let us define 

oo 

v(k,t) = J2 e~ ikn u n (t), (4) 

n=— oo 

where u n (t) is related to v(k,t) as follows 

TV 

u n {t) = — I dke ikn v(k,t), (5) 

Z7T J 

— n 

and k can be considered as a wave number. 
Then the second term in Eq.(^) reads 

1 1 7 

- V n - m (u n -u m ) 2 = — / dkV(k)\v(k,t)\ 2 , (6) 

where function V(k) has been defined by 

V(k) = V(0)-V(k), (7) 



and V(k) is 



V{k)= £ e~ ikn V n . (8) 



n=— oo 
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In the long wave limit when the wavelength exceeds the intersite scale we 
may consider v(k, t) as a fcth Fourier component of function u(x, t), u n (t) - — > 

u(x, t). That is the functions v(k, t) and u(x, t) are related each other by the 
Fourier transform 



1 °° 

u (x,t) = — J dke lkx v(k,t), v(k,t)= I dxe- ikx u(x,t). (9) 



Therefore, from Eq.fJH} we have 



1 - - , 1 



^ 51 v n-m(u n - u m ) 2 = - J dx J dyu(x , t)W (x - y)u(y , t) , (10) 



^ n— — x hi — — x 



with the kernel 



—00 —00 



00 

1 



W(x) = - / dke lkx V(k). (11) 

7T 



— 00 



that comes from the long-range interaction matrix V n _ m and it possesses the 
properties: 

(i) W(x) is an even function; 

00 

(ii) / dxW(x) = 0. 

—00 

One can express Eq.([10|) in a non-local kinematic form 



00 00 

1 ^ , . ,0 1 



- J n-m(u n - u m ) 2 = - f dx j dy(d x u(x,t))K(x-y)d y u(y,t), (12) 



n,m 

_/ —00 — 00 



where notation d x = d/dx has been introduced and the relationship between 
kernel K(x) and function V(k) (see Eq.flJJ)) is given by 

K(x) = - J dke lkx ^-. (13) 

—00 

In other words, the kernels W(x) and K(x) are related to each other as 
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W(x) = -d 2 x K(x). (14) 
Thus, in continuum limit the equation of motion becomes 

7 dU(u) 
Mu{x, t)- J dyd x K{x - y)d y u(y, t) + — A > = 0, (15) 

— oo ' 

or in a symbolic form 

Mii(x, t) - d x (Kd x u) + Q^jj = 0, (16) 

and the operator K, as it follows from Eq. (JTHJ), can be considered as operator 
of multiplication in wave number (momentum) space 

Ku(x,t)(k) = ^-v(k,t). (17) 

Equation (|15j) is integro-differential equation of motion. The integral part 
comes from the long-range interaction term Eq. (jl2|) . To get the differential 
equation of motion we use the properties of function V(k) at k — > 0, which can 
be obtained from the asymptotics of polylogarithm (we provide the properties 
of the polylogarithm function and the expression for V(k) in the Appendices), 

V(*)~— ^^Ikr 1 = D a \k\'~\ 2<s<3, (18) 



T(s)sin( , 



V n k 



2 



V(Jfe) ^lnfc 2 , 8 = 3, (19) 

V(k)r,Y^i^le, s>3, (20) 

where T(s) is T-function, ((s) is Riemann zeta function and coefficient D s is 
defined by 

D a = ^W- (21) 

r( s ) sin (E(^i) 1 ; 

It is seen from Eq. ()18|) that fractional powers of k occurs for the interactions 
with 2 < s < 3 only 

Going back to the coordinate space yields for fractional powers of \k\ 
fractional Riesz derivative of order a jl^j 
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1 / d a d a \ 

8y(x, t) = , — — + — — u(x, t), l<a<2, (22) 

x v ; 2cos(7ra/2) \dx a d(-x) a J v ' ' - , \ j 

where d a /d(±x) a are the Riemann-Liouville derivatives 

d a 1 d n J dyf(y) 

~J\ X ) 



dx a v ' T(n-a)dx n J (x - y)°- n+1 



oo 



d a _ f{x) _ I d n f dyf(y) 



d(-x) a T(n-a)dx n J x (y - x) a ~ n+1 ' 

with n — 1 < a < n and integer n. Substitution of Eq. (fT£j) into Eq. (fTH|) yields 
a fractional nonlinear wave equation 

Mu(x, t) - D s d s ~ l u{x, t) + = 0, (23) 

where the Riesz fractional derivative has been transformed into 



<9*-yx,t) = / ^^^^-^(fc^), (24) 

27T 7 

— oo 

and coefficient D s is defined by Eq. (j21)) . 

When the nonlinear external potential U(u) = Uq(1 — cosu), we get the 
fractional sine-Gordon equation 



Mii(x, t) - D s d s - l u(x, t) + U Q sin u(x, t) = 0. (25) 

To get differences for the phonon spectrum of fractional and standard sine- 
Gordon equations let's consider a linear equation of the form 

il(x, t) - c 2 d 2 xx u{x, t) + m 2 u(x, t) = 0, (26) 

where cq is velocity of wave and m is the particle mass. Plane wave solution 
of this equation has the form 



■ft 



(x,t) = Ae i{Mx \ (27) 
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with constant A, frequency oj, and wave number k. Then the dispersion law 
is 



u(k) = ±^Jm 2 + elk 2 . (28) 
For the fractional linear wave equation 

ii(x, t) - b s -*cld s - x u{x, t) + m 2 u{x, t) = 0, (29) 

where d s x ~ l is defined by Eq.(J2U), b is the scale constant attributed to the 
fractional wave equation, c is wave velocity, the dispersion law is 



fph k 



u(k) = ±^m 2 + b s - 3 c 2 \k\ s ~\ 2 < s < 3. (30) 
For fairly small values of m and not too small k we obtain 

co(k) ^M s - 3 >/ 2 c |£;| (s - 1)/2 , (31) 
which results in phase v p h and group v g velocities equal 

^1 = b^ 2 c /\k\^l\ 2<s<3, (32) 
k 

v 9 = ^ = ^- 3)/ Vl^l (3 - s)/2 = ^v p ,, 2< S <3. (33) 

These expressions tend to infinity for k — > 0, and we arrive at new physical 
properties of the lattice of particles just because of fractality of the long- 
range interaction. For small enough k the mass m can not be neglected if 
m^O. Several other examples of equations can be easily obtained using the 
general scheme of Eq. ()23|) 

It follows from Eq.()18|) that in the case s = 2 the function V(k) takes the 
form 

V(jfe) = 7rVb|A;|. (34) 

Substitution of this expression into Eqs.(|13|) and (|16j) yields a nonlinear wave- 
Hilbert equation 

Mu(x, t) - 7rV H {d x u(x, t)} + = 0, (35) 

ou[x, t) 

where H. is the Hilbert integral transform defined by 
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H{<P(x)} = P J dy^-, (36) 

— oo 

and P stands for the Cauchy principal value of the integral. 

When U(u) = U (l — cosw) we find the sine-Hilbert equation [31], [32] 

Mu(x,t) — 7rV H {d x u(x, t)} + U smu(x,t) = 0. (37) 

Examples of physical problems with nonlocal interaction like V(k) given 
by Eq.fjHD) can be found in nonlocal Josephson electrodynamics [3^1 where 
Eq.()37|) is one of the basic model equations. Let's note that some exact 
solutions of this equation have been found and their stability under weak 
perturbations has been analyzed [33]-|3T]. 

Finally, when s > 3 and U(u) = Uq(1 — cosw) we get the sine-Gordon 
equation [3J, 

Mu(x, t) - V °^~ 2 ^ dl x u(x, t) + U sin u(x, t) = 0, (38) 
where y ^tr' 2 ^ can be interpreted as the wave velocity. 



3 Quantum lattice dynamics 

3.1 Lattice Hamiltonian with long-range interaction 

To model a ID quantum lattice dynamics we follow [5], [S] and consider 
a linear, rigid arrangement of sites with one molecule at each lattice site. 
To describe creation (annihilation) of molecular excitation or, for simplicity 
exciton, at the site n we introduce exciton creation 6+ and annihilation b n 
operators. Operators 6+ and b n satisfy the commutation relations [b^,b m ] = 
$n,mi [b n , b m ] = 0, [&+, £>+] = that is 6+ and b n are the Bose operators. 

With the help of the operators 5+ and b n the well-known Hamiltonian 
operator of excitons is expressed as [Sj, [E], |3~Hj . 



H D = £ X] Kbn-J Yl {b+b n+ i+b+b n -i) = e ^ b+b n - ^ J^mb+bm 

(39) 



n,m=— oo 
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where e is constant exciton energy, J is interaction constant and the exci- 
tation transfer matrix element J n>m between sites n and m can be written 

as 

+ 5 

(n— 1), m ), (40) 

the Kroneker symbols 5 m ,n±i mean that only the nearest-neighbor interaction 
had been considered. The interaction term J n ,mbtb m in Eq.()39j) is responsible 
for exciton transfer from site n to the nearest neighbor sites nil. 

To go beyond the nearest-neighbor interaction we introduce the long- 
range excitation transfer matrix which describes power-law interaction 
between sites n and m (see Eq.©) 

J L n % ee V n . m = n ± m. (41) 

Then we can rewrite the Hamiltonian operator H^r with long-range excita- 
tion transfer, 

Hlr = e £ b + n b n - £ ( 42 ) 

Follow Davydov's ansatz [5], [U] we define the eigenstate \4>{t) > of a 
quantum exciton as a superposition: 

|0(*)>=EV'»(*)6J|O>, (43) 

n 

where 1 > is the vacuum state of exciton system and ip n {t) is the exciton wave 
function. In the \4>(t) > basis the Hamiltonian Hlr becomes the Hamiltonian 
function H LR (ip, ip*) 

h lr (w) =< m\H LR m) >=ej2r n m n (t) - e k^cw^w. 

n n,m 

(44) 

where V n - m is given by Eq.(jHJ). Considering nonlinearity in addition to the 
long-range interaction, we add a nonlinear term C/(|^|) to the right hand of 
Eq . (|44j) . and, thus, we obtain for the Hamiltonian H(tp,ifj*) of ID lattice 

%f) = £ Ecwtw- e Vn- m r n m m (t)+Y, u (\ipn\)- (45) 

n n,m n 
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The nonlinear Schrodinger equation in discrete space ifidifj n (t) / dt = SH(ip)/S^ 
generated by the above Hamiltonian has the form 

0*W!> =*.(,)- ■£ + (46) 

at v d%p n 

where % is the Planck constant. 

To get continuum space equation let us follow the same approach as in 
Sec. 2 and introduce wave function (p(k,t), 

oo 

<p{k,t) = £ e" ifc > n (t), (47) 

n=— oo 

where ip n (t) is related to <f(k,t) as 

1>n(t) = 7r I dke ik Mk,t), (48) 

— IT 

and k can be considered as a wave number. Substitution of Eq. ()48)1 into 
Eq.(gU) yields 

1 n 

E lr = — \ dk(u + V(k))\ip(k,t)\ 2 , (49) 

with the function V(k) given by Eq.© and the energy parameter u which is 

u = e-V{0), (50) 

oo 

here ^(0) = Vn (see the definition given by Eq.(JHJ)). 

n=— oo 

At this point we may consider tp(k,t) as the wave number representation 
of continuous space wave function ip(x,t). That is the functions (p(k,t) and 
ip{x,t) are related to each other by the Fourier transform 



1 x 



${x,t) = — / dke lkx <p(k,t), (p{k,t)= dxe~ lkx ilj(x : t). (51) 

—oo 

Therefore, we have for Hlr in coordinate space 
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oo oo 

,2 1 



H LR {^)=u I dx\^(x,t)\ 2 + - / dx I dy^*(x,t)W(x-y)^(y,t), (52) 



-oo — oo 



where the kernel W(x) is defined by Eq.(jllj). To get kinematic form of the 
second term in the right side of Eq. (|52|l we follow Eqs.lfT^j) and (jHj) and find, 



oo oo 



H LR (i))=uj I dx\^(x,t)\ 2 + - I dx I dy(d x ^*(x,t))K(x-y)d y ^(y,t), 



-oo — oo 



(53) 

where the kernel K(x) is defined by Eq.(|T3*|). 

Further, by adding a nonlinear term to right hand of (|53|1 we obtain the 
integro-differential (non local) nonlinear Schrodinger type equation which 
can written in the following compact way 

a ^ = ^ +M + %ML y (54, 
where the operator K acts on wave function ip(x,t) as 

K1>(x, t) = ^J dke ikx ^-v(k, t), (55) 

— oo 

or in wave number (momentum) representation 

K^(x,t)(k) = ^-p(k,t). (56) 

Using the properties of function V(k) in the limit k — > (see Egs- lfTSjl - lj^Ojl ) 
we can get different special forms of general Eq. (j54|) . 

3.2 Fractional nonlinear Schrodinger equation 

Substitution Eq.flTHJ) into Eqs.(|5^|l and (|53j) yields a new nonlinear fractional 
differential equation 



in ^EA = ^ t) _ Dsd s-i^ t) + EML } 2<s<3, (57) 
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where we use the definition of the Riesz fractional derivative j2H] > [2*B] 



3C 



fi£-ty(:M) = -— / dke^lkl 3 - 1 ^,^. (58) 

— oo 

In Eq.(|57j) coefficient D s is given by Eq. (|21|) and uj is the constant defined 
by Eq.([50j). For the case when 

U(\i;\) = U(xM 2 , (59) 
Eq. (JH7|) is reduced to the linear fractional Schrodinger equation [32] , [EE] , j3D] 

.^^f) = _ Dsd s-i^ t) + u( x)< f,( Xi ^ (60 ) 

where the wave function 4>{x, t) is related to the wave function ip(x, t) by 

<p(x,t) = exp{i—}ip(x,t). (61) 

Generalization of Eq. (j60|) for 3D space and the applications of the ID and 
3D linear fractional Schrodinger equation to quantum mechanical problems 
have been developed in jUJJ. In these papers, a particle in infinite potential 
well, fractional oscillator, and fractional Bohr atom have been studied and 
the energy spectra for these three quantum mechanical problems have been 
obtained using Eq.([60|) and its 3D generalization. New physical issues follow- 
ing from quantum mechanical applications of Eq. (j60p have been discussed in 

EH- 

For nonlinearity of the form 

UM) = ^H 4 , (62) 

with a constant Uq, Eq. (jBT|) is reduced to the fractional nonlinear Schrodinger 
equation 



m d^t) = _ Dsd s-i (f) ^ t) + t)] 2^ t)? (63) 

where <p(x, t) is related to the wave function ip(x, t) by means of Eq. (J61|) . The 
equation and its 3D generalization were proposed in jH] to study wave 
propagation or kinetics in a nonlinear media with fractal properties (see also 
[4*2*] . |43j). Following from Eq. (p)3*j) a fractional generalization of the nonlinear 
Ginzburg-Landau equation has been developed and square integrability of its 
solution has been established as well in |4*T] . 
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3.3 Nonlinear Hilbert-Schrodinger equation 

For the case when s = 2 substitution Eq.(fT%|) into Eqs. ljBTj) and (JS3j) yields a 
new nonlinear integro-differential equation 

irfWzA = ^ t) _ m{x ^ t)} + (64) 

where Tt is defined in Eq. (jSEJ) 

When Z/dV'l) is given by Eq. (|59jl we find from Eq. (|67|) the linear Hilbert- 
Schrodinger equation 

if d^t)_ = _ W ^ H t)} + u{x)(j){x ^ f y (65) 

When tZ(lV'l) is given by Eq.([62p we obtain from Eq.()67|) the nonlinear 
Hilbert-Schrodinger equation 

iU 8 ^' = -^V H{d x 4>{x,t)} + U Q \4>{x,t)\ 2 (j>{x,t). (66) 

Finally, we note that for s > 3 it follows from Eq.()2L)|) that Eq.(j57|) takes 
the form of the nonlinear Schrodinger-like equation with extra term LVip(x, t) 
originated because of the energy parameter uj 



here is the Riemann zeta function. 



4 Quantum exciton propagator 

To get insight on impact of long-range interaction on ID quantum dynamics 
let us focus on the discrete linear problem associated with Eq. (j46|) 

^ = «- E V n ^ m (t), (68) 
with long-range potential V n - m given by Eq.(j2*|). 
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Suppose that we know the solution ip n >(t') of Eq.(|55jl at some time instant 
t' at the site n . Then solution ip n (t) at later time t, (t > t'), and site n will 
be 

1p n (t)=Y, G n,nW)lPn'(t'), (69) 
n' 

where G n>n i(t\lf) is quantum exciton propagator, that is probability of exciton 
transition from site n' at the moment t' to site n at the moment t. 

It follows from Eq.(|B5|l and Eq. fJBTJj) that G n . n '{t\t') is governed by the 
equation 



ih dG n At\t>) =eGnAAf) _ ^ V n „ m G m , n ,(t\t'), t>t\ (70) 

m 

with the initial condition G n ^ n i(t\t) = 5„ )n ', here 5 n>n ' is the Kronecker symbol. 
Let us put for simplicity n' = and if = and introduce a notation 

G n (t) = G n>0 (t\0). (71) 

It yields 

l% dGM = eGn{ t) - £ V n - m G n (t), t>0, (72) 
ot m 

with the initial condition 



G n (t)\ t=0 = 5 nfi . (73) 

To get the continuum version of Eq. (J72)) we apply transformations similar 
to Eq.(jlJ) and Eq.©. That is we define G(k,t) as 

oo 

G(k,t)= J2 e- lkn G n (t), (74) 

n=— oo 
1 ^ 

Gn{t) = — f dke lkn G(k,t), (75) 

271 J 

— 7V 

and G(k,t) satisfies the equation 



16 



ih^^- = {u + V{k))G{k,t), t>0, (76) 

where V(k) is the same as in Eq.((7j) and the energy parameter uj is defined 
by Eq. (JHO|) . The initial condition Eq. (J73Jl now becomes 

G(k,t = 0) = I. (77) 

For further consideration it is convenient to introduce the quantum propa- 
gator g(k,t) which is 

g(k,t) = exp{i^}G(k,t). (78) 
Thus, G(k,t) can be expressed in terms of g(k,t) as follow 

G(k,t) = exp{-i^}g(k,t). (79) 
It follows from Eq.([7fi|) that the propagator g(k,t) is governed by 

z^j^ = V(%(M), *>0, (80) 
with the initial condition 

g(k,t = 0) = 1. (81) 
The solution of the problem Eqs. (j80j) and (JHTj) is 

g(k,t) = exp(i-^-). (82) 
By substituting this solution into Eqs.(|79^1 and (|75j) we find 

G n {t) = j^J dkexp(tkn - i ( u + VW\ (83) 
or, after extraction of the energy-time factor exp{— i^}, 

G n (t) = exp{-Ag n (t), (84) 
where the lattice propagator g n (t) has been introduced by 
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9n(t) = t-^t J dkexp(ikn - i —^- ). (85) 

— 7T 

Generalization to ID lattice quantum exciton propagator g n y(t\t') which 
describes transition from site n' at the time instant t' to site n at the time 
instant t, is obvious 

g n ,n>W) = J^J dkexp{tk{n - n') - j ^^"^ }. (86) 

— 7T 

This is the quantum exciton propagator which describes ID transport discrete 
in space and continuous in time. Let us note that g n , n '(t\t') satisfies the 
following conditions: 

1. Normalization condition 

oo 

£ g n>n ,{t\t') = 1, (87) 

n=— oo 

2. Consistency condition 

oo 

9n,n'(ti\t 2 ) = 9n,m(tl\t') ' g m ,n' (t'\t 2 ) ■ (88) 

m=— oo 

The last condition means that exciton propagator g n ,n'(t\t') can be con- 
sidered as a transition quantum amplitude, and for propagations occurring 
in succession in time transition amplitudes are multiplied. 

Further we will study behavior of g„(t) at large \n\, when the main con- 
tribution to the integral Eq. (J85|) comes from small k. Therefore, we expand 
integral over k from — oo up to oo 

k< s h. 



with 



and 



I r t k u(s 't 

9n(t) = TT / dkexp{ikn - i^-r—), (89) 
2tt J n 



( x r 2 , for s > 3, 
" (a) = { a -l, for2<,<3, (90) 
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^C(s-2) for s > 3 
ls 1 D s , for2<s<3, 1 1 

where D s is given by Eq. (j21|) . 

Asymptotic behavior of the quantum ID propagator g n (t) at large |n| 
depends on the parameter s. Indeed, when s > 3 Eq. (|H9*j) goes to 



% I 1 2 

B (f) = (ft/27rzV C(s - 2)t) 1/2 exp{- 2 ^ c ^_ 2)t } ; S > 3, (92) 

if we take into account Eqs. (|9"Uj) and (jUT|) . 

When 2 < s < 3 the integral in right-hand side of Eq. (|H9~|) can be ex- 
pressed in terms of the Fox's if-function [44] and we have 



(l,l/( S -l),(i,l/2)\ 
(1,1), (1,V2) J ' 

(93) 

where if 2 ' 2 is the Fox's function (for definition see, for example [13], [45J). 
From other hand, in the long- wave limit for 2 < s < 3 the integral in Eq. (j89j) 
can be estimated as 



9n{T) 



\n\{s — 1 



-H, 



i,i 

2.2 



n 



i/(»-i) 



iD, 



n 



T 



9n (t) ~ -r(a) sin(^-^) • — , 2 < s < 3. (94) 

Thus, the long- wave asymptotics at large \n\ of the quantum exciton propa- 
gator g n (t) exhibits the power-law behavior for 2 < s < 3. Transition from 
Gaussian-like behavior Eq. (J92j) to power-law decay Eq . (J94j) is due to long- 
range interaction (second term in the right hand of Eqs. (j68|) . (J72|) ). This 
transition can be interpreted as phase transition. In fact, for s > 3 propa- 
gator g n (t) given by Eq. ()92|) has correlation length (or characteristic scale) 
which can be describe by 

f 2V qs-2)t \ 1/2 

An ^{ — h — J ' 

as far as in all above considerations we put the lattice constant a equal 1 for 
simplicity. As it follows from Eq . (J94j) . the correlation length is infinite, that 
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is the correlation length doesn't exist for long-wave excitons in ID lattice 
with 2 < s < 3 power-law interaction. 

Finally, let us note that the above discussed ID quantum lattice dynamics 
for imaginary time is similar to the random walk model, that is similar to 
CTRW 24 . Indeed, if we put it — > r, Eq.fl8H|) is transformed to 

Pn,n'{r\r') = T^y / dkexp{ik(n - n') - YBi^Zll }, (95) 

— 7T 

where P n ,n'{ T \ T ') is transition probability, i.e. the probability to walk from 
site n to site n' at the time interval r — r'. It is easy to see that P n ,n'( r \ T ') 
is normalized 

00 

E PnAr\r') = l, (96) 

n=— 00 

and satisfies the discrete version of the Kolmogorov-Chapmen equation 

00 

iV(nk 2 )= E ^n, m (ri|r / )- J P mX (r'|r 2 ). (97) 

m=— 00 

That is we came to a continuous in time and discrete in space (ID lattice) 
random walk model. The obtained random walk model exhibits a phase 
transition from the Brownian random walk (s > 3) with finite correlation 
length to the symmetric a-stable (with a = s — 1, 2 < s < 3) or the Levy 
random process with an infinite correlation length. 

In the case of continuum space the integral over k Eq. (|95|) can be ex- 
panded from —00 up to 00 and we obtain transition probability distribution 
function P(x, t\x', t') 

P(x,t\x',t') = (98) 

— J dke*pUk(x-x) j, 

—00 

where v(s) and j s are defined by Eq. (j9()j) and Eq.(|91jl reciprocally. This 
probability distribution function exhibits transition from Gaussian-like be- 
havior at s > 3 to power-law decay at 2 < s < 3 because of the power-law 
long-range interaction. 
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5 Conclusions 



We have presented analytical developments on classical and quantum ID lat- 
tice dynamics with fractional power-law long-range interaction. The classical 
case has been treated in the framework of the well-known Frenkel-Kontorova 
chain model. Quantum lattice dynamics is considered following Davydov's 
approach to dynamics of molecular excitons. It has been shown that the long- 
range power-law interaction leads, in general, to a nonlocal integral term in 
the equation of motion if we go from discrete to continues space. In some 
particular cases with non-integer power s for power-law interaction the inte- 
gral term can be expressed through the fractional order derivative. In other 
words nonlocality from the long-range interaction reveals the dynamics in 
the form of space derivatives of fractional order. 

As useful examples of classical lattice dynamics in the continuum media 
limit we were able to derive in a unified form fractional sine-Gordon and 
wave-Hilbert nonlinear equations. 

In the quantum case we obtained linear a fractional Schrodinger equa- 
tion, fractional nonlinear Schrodinger equation, linear Hilbert-Schrodinger 
equation, and nonlinear Hilbert-Schrodinger equation. 

This applied approach permits one to see, in an explicit form, the in- 
terplay between second order and fractional order space derivatives both in 
classical and quantum dynamics. 

We have observed that there exists critical value for the power of the long- 
range potential such that separates the cases of fractional equations from the 
well-known nonlinear dynamical equations with integer order derivatives. 

The long-range interaction impact on quantum lattice propagators has 
been studied. It was shown that the exciton propagator exhibits a transition 
for the same critical value of the power-law exponent from the well-known 
Gaussian-like behavior to a power-law space decay. This transition can be 
treated as a phase transition accompanying by an infinite growth of the 
exciton correlation length. 

It has been shown that in the imaginary time domain ID quantum lat- 
tice dynamics with long-range interaction can be considered as random walk. 
The corresponding equations (Eq.(|95|l and Eq.(|9H|) ) for the transition prob- 
ability of this random walk have been obtained. The random walk model 
also exhibits a phase transition from the Brownian random walk (s > 3) 
with finite correlation length to the Levy random walk (2 < s < 3) with an 
infinite correlation length. 
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Finally, let us note that this approach can be generalized to 2D and 3D 
classical and quantum lattice dynamics. 
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7 Appendix 

7.1 Poly logarithm as a power series 

The polylogarithm Li s (z) is defined by 



oo ~n 7 t f s ~ 1 



n=l 



here s is real parameter and argument z is the complex argument. It is easy 
to see that for z — 1 the polylogarithm Li s (l) reduces to the well-known 
Riemann zeta function ((s), 

OO 1 

Lis(l) = C{s) = £-■ (100) 

n=l U 

It seems that the quantum statistical mechanics is the best known field where 
the polylogarithm arises in natural way. Indeed, the integral of the Bose- 
Einstein distribution is expressed in terms of a polylogarithm as, 

OO _^ 

Li 'W = fw/ dt i7^T' (101) 

v ' ' 

In statistical mechanics handbooks the above integral is referred as a Bose 
integral or a Bose-Einstein integral (see, for example |IHj)- 

Next, the integral of the Fermi-Dirac distribution is also expressed in term 
of a polylogarithm, 
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u -^ = ~W) J dt ^NTT- (102) 

v ' ' 

This integral is referred as a Fermi-Dirac integral [4T?j . 

It is convenient to write the argument as e^, that is consider Li s (e M ). 
Then one can use the Hankel contour integral for the polylogarithm [571] . 

Li^) = -^^/^^4, (103) 

H 

where H represent the Hankel contour and s ^ 1,2,3,..., and pole of the 
integrand does not lie on non-negative real axis at the t = fi. The integrand 
has a cut along the real axis from zero to infinity, with the real axis being 
on the lower half of the sheet (Imt < 0). For the case where fi is real and 
non- negative, contribution of the pole has to be count, 

Li s (e^) = <j> dt^-j-jf—- + 2niR, (104) 

Z7T J e l v — 1 
H 

where R is the residue of the pole 

rq 

2tt V ' 

To get power series (about fi — 0) representation of the polylogarithm 
Li s (e~ tl ) we apply the Mellin transform following Ref. 



M s (r) = J duu r - l U s (e- u ) = — J du J dt^—j. (106) 

Then we observe that changing of variables t = xy, u = x(l — y) separates 
the integrals and we obtain 



M s (r) = -L| dyy'-^l-y)- 1 ] dx^-^ = r(r)Li s+r (l) = r(r)C(* + r). 
^ ' o o 

(107) 

Now through the inverse Mellin transform we have 
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c+ioo 

Li s (e^) = — / dr/i"T(r)C(s + r), (108) 
2m J 

c—ioo 

here c is a constant to the right of the poles of the integrand. The path 
of integration may be converted into a closed contour, and the poles of the 
integrand are those of the gamma function T(r) at r = — I with residues 
(— 1) //! (I = 0,-1,-2,...), and of the Riemann zeta function q(s + r) at 
r = 1 — s with residue +1. Summing the residues yields, for < 2n and 
s ^1,2,3,.... 

Li s (e^) = r(l - s)^)- 1 + £ ii—J^y. (109) 

This equation gives us is get power series (about fi = 0) representation of 
the polylogarithm Li s (e~ M ). 

Let's note that if the parameter s is a positive integer n, both the I — n— 1 
term and the gamma function T(l — n) become infinite, although their sum 
does not. For integer I > we have 

^i-H (110) 

and for / = 0, 

lim [((s) + r(l - S )(-M) S_1 ] = - In(-Ai). (HI) 

s— >1 

The power series representation Eq. (jl09j) will be used to get Eqs. (lll4|) - ()116j) . 



lim 



7.2 Properties of function V(k) 

It is easy to see that V(k) given by Eq. (JSJ) with V n defined by Eq.Q can be 
expressed in terms of the polylogarithm Li s (z) 

V(k) = 2V f) = 2V ■ Re{Li s (e- lfe )}. (112) 

Further, the function V(k) defined by Eq.© becomes 
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V(k) = 2V ((s)Re(l- U jf* k) ). (113) 

Taking into account power series representation given by Eq. (jl(J9j) we find 
from Eq. ([113)1 at the limit k — > that 

V(k)~ - ; V %- lh \kr\ 2<s<3, (114) 
T(s) sm( v 2 ' ) 

V(fc) ~ -^lnfc 2 , s = 3, (115) 

V(fc)~ y ° C(a ~ 2) fc 2 , S >3. (116) 
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